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The electronic structures of Ba(Fe,Ru)2As2 and Sr(Fe,Ir)2As2 are investigated using density func- 
tional calculations. We find that these systems behave as coherent alloys from the electronic struc- 
ture point of view. In particular, the isoelectronic substitution of Fe by Ru does not provide doping, 
but rather suppresses the spin density wave characteristic of the pure Fe compound by a reduction 
in the Stoner enhancement and an increase in the band width due hybridization involving Ru. The 
electronic structure near the Fermi level otherwise remains quite similar to that of BaFe2As2. The 
behavior of the Ir alloy is similar, except that in this case there is additional electron doping. 

PACS numbers: 74.25. Jb,71.20.Be 



I. INTRODUCTION 

The discovery of unconventional superconductiv- 
ity in proximity to magnetism for layered Fe-based 
materials 1 ^^^ has led to considerable interest in both 
establishing the interplay between magnetic order and 
superconducting state, and searching for effective ways 
of tuning them. The phase diagrams typically show a 
spin density wave (SDW) that competes with supercon- 
ductivity, i.e. superconductivity generally appears when 
the SDW is suppressed either by doping with electrons 
or holes, which reduces the nesting by making the elec- 
tron and hole Fermi surfaces mismatched or by pressure — 
which increases hybridization and broadens the bands 
again working against nesting. 

A remarkable feature of these compounds is that, in 
contrast to the cuprates, they may be doped into a su- 
perconducting state by substitutions on the Fe-site, e.g. 
in Ba(Fe,Co) 2 As 2 and Ba(Fe,Ni) 2 As 2 ££ In these alloys 
the electronic structure remains very similar to that of 
the parent Ba(Sr)Fe 2 As 2 , but the Fermi level is shifted 
upwards corresponding to the increased electron count. 
Thus one difference from cuprates is that substitution 
of Co and Ni lead to the formation of a coherent alloy 
electronic structure rather than the introduction of lo- 
calized states associated with those ions. Recently, it 
has been shown that in addition to the 3g? dopants, Co 
and Ni, superconductivity can be induced by alloying 
with some 4c? and 5c? elements, including Ru, Rh, Ir and 
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Transition elements in the 4c? and 5c? rows differ from 
their 3c? counterparts in several respects. Since Ad and 
5c? orbitals are much more extended than the 3d orbitals 
there is a greater tendency towards covalency both with 
ligands (e.g. As) and also in stronger metal - metal d 
bonds. For example, mid-5c? transition elements have 
some of the highest melting points of any material (the 
melting points of Ir and Ru are 2739 K and 2607 K, re- 
spectively, as compared to 1811 K for Fe and 3695 K 
for W), and compounds of these elements are often ex- 
tremely hard. Thus one may expect broader bands and 
more hybridization with As in the alloys with these el- 
ements. Secondly, again because of the more extended 



4c? and 5c? orbitals, the Hund's coupling on these atoms 
is weaker than on 3c? atoms, which works against mag- 
netism and is reflected in lower values of the Stoner pa- 
rameter for 4c? and 5c? elements and compounds. Finally, 
the larger orbitals lead to a tendency for higher valence 
states in the 4c? and 5c? series, so that Ru 4+ and Ru 5 + 
compounds are more stable and more common than the 
corresponding Fe compounds. 

Returning to superconductivity in alloys of BaFe 2 As 2 
with 4c? and 5c? elements, the case of Ru is particularly 
interesting because Ru is isoelectronic with Fe and there- 
fore it may or may not serve as a dopant. Transition tem- 
peratures up to T c ~ 21 K may obtained with substantial 
nominal Ru content in Ba(Sr)Fe 2 _ x Ru x As 2 , x nom ~0.7 
although Ba(Sr)Ru 2 As 2 shows neither magnetic order 
nor superconductivity^ One possibility is that Ru serves 
as a dopant, supplying carriers to the Fe planesi 10 ' 12 For 
example, Ru might occur as Ru 4+ , in which case empty 
localized d states associated with Ru atoms would oc- 
cur above the Fermi energy, while the Fe derived valence 
bands would show higher filling reflecting electron dop- 
ing by two carriers per Ru. Another possibility is that 
alloy shows a coherent electronic structure that is dis- 
torted from the pure Fe compound, but does not show 
additional states, similar to the Co and Ni doped materi- 
als. In this case, the effect of Ru could be similar to that 
of pressure, broadening the bands, increasing hybridiza- 
tion and/or lowering the effective Hund's coupling and 
thereby destroying the SDW in favor of superconductiv- 
ity. 

Here, we report density functional studies of the elec- 
tronic structure of BaFe 2 _ a; Ru a ;As 2 . It is found that the 
substitution of Fe with Ru gives a quite similar electronic 
structure to that of BaFe 2 As 2 near the Fermi level Ep. 
Importantly, this substitution does not result in doping 
since we find no additional flat bands reflecting Ru states 
in supercell calculations and consequently the exact Lut- 
tinger's theorem compensation of electron and hole sur- 
faces is maintained in the same way as in the pure com- 
pound without the SDW. Thus, Ru substitution does not 
introduce additional electrons to this system. Instead 
the SDW magnetic order is suppressed by the decreasing 
Stoner term and increased hybridization. We also per- 
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formed the calculation for 5d transition metal, Ir alloyed 
SrFei_ x Ir a ,As2, which was reported to superconduct with 
T c ~ 22 K. 1 — We also find an electronic structure charac- 
teristic of a coherent alloy for this non-isoelectronic, Co 
column transition metal. We find electron doping as ex- 
pected with additional one electron per Ir, which cooper- 
ates with the reduction of Hund's coupling to destabilize 
the SDW order. 



II. FIRST PRINCIPLES CALCULATIONS 

Our first principles electronic structure calculations 
were performed within the generalized gradient ap- 
proximation (GGA) of Perdew, Burke, and Ernzerhof 
(PBE), 17 using the general potential linearized aug- 
mented planewave (LAPW) method, with the augmented 
planewave plus local orbital implementation i 18 ' 19 i 20 
LAPW spheres of radii 2.3ao for Ba, Ir and 2.1ao for 
Fe, As, Ru were employed. We used well converged basis 
sets determined by -R m i n -K max =8.5, where i? m i n is the ra- 
dius of the smallest LAPW sphere and K max corresponds 
to the planewave cutoff for the interstitial region. We in- 
cluded relativity at the scalar relativistic level. Local 
orbitals were included to accurately treat the semi-core 
states. The Brillouin zone sampling for self-consistent 
calculations was done using the special k-point method, 
with a 24x24x24 grid for body-centered ThCr2Si2 struc- 
ture and a 21x21x9 grid for the quadrupled tetragonal 
supercells with partial Ru and Ir substitution (see be- 
low). 
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FIG. 1: Calculated electronic band structure of (a) BaRu2As2 
with the relaxed zas, comparing with that of (b) BaFe2As2. 
The bands are plotted along high-symmetry directions in the 
body-centered tetragonal Brillouin zone, and X point corre- 
sponds to the zone corner M point in the tetragonal Brillouin 
zone. For BaFe2As2, we used the experimental lattice pa- 
rameter a = 3.9625 A and c = 13.0168 A (Ref. H) and 
non-spin-polarized GGA optimized As height zas = 0.344. 



III. RUTHENIUM SUBSTITUTION 

We first calculated the electronic structure of com- 
pletely Ru substituted non-superconducting BaRu2As2 
to check whether additional bands or significant changes 
in the band shapes are introduced by Ru. The experi- 
mental tetragonal lattice constant a = 4.152 A and c = 
12.238 A 21 were employed and the internal As coordi- 
nate zas were optimized by total energy minimization as 
0.3528. For this material, we checked for but did not find 
magnetic order of cither ferromagnetic or antiferromag- 
netic character, consistent with recent experiments. 16 
Fig. [T] shows the calculated band structure, in compari- 
son with that of non-spin-polarized BaFe2As2 as obtained 
with the same approach. As may be seen, the general fea- 
tures of band structures around Ep for two compounds 
show quite close similarity and no new bands appear in 
BaRu2As2- In particular, there are similar compensat- 
ing electron Fermi surface sections at the zone corner 
and hole sections at the zone center (r-Z direction). Es- 
pecially at the zone corner, there is a very strong sim- 
ilarity of the electron cylinders in the two compounds. 
However, for BaRu2As2 we do find a somewhat different 
structure of the hole Fermi surfaces, which are centered 
near Z as compared to BaFe2As2 where hole bands exist 



along the T-Z direction. This yields a closed hole Fermi 
surface in the k z direction in BaRu2As2- This more 3D 
shape is suggestive of more Ru - As hybridization, which 
is also seen in the projections of the electronic density 
of states (DOS). In any case, the position of Ep with 
respect to the compensating position between the hole 
and electron band edges is maintained similar to non- 
spin-polarized BaFe2As2. This is different from the non- 
isoelectronic compounds, BaCo2As2 2 ^ and BaNi2As2r^ 
which also show rather similar electronic structures to 
the Fe compound, but with the Ep shifted to higher en- 
ergy corresonding to the higher electron count. This shift 
of the Fermi energy results in dramatically different phys- 
ical properties. 

We used supercell calculations to explicitly study the 
effects of partial Ru and Ir substitution. For this purpose 
we used a \Z2xV2 doubling in plane, and also doubled the 
cell along the c-axis by going to the conventional tetrag- 
onal cell as opposed to the body centered cell. This leads 
to a quadrupling of the ThCr 2 Si2 structure cell. We then 
replaced one Fe by Ru or Ir in each plane. This cor- 
responds to a composition of BaFe1.5Ruo.5As2 (nominal 
Ba4Fe 6 Ru 2 As 8 ) or SrFe1.5Iro.5As2. For Ba4Fe 6 Ru 2 As8 
we did calculations for three different arrangements of the 
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FIG. 2: (Color online) Comparison of the calculated electronic 
DOS for the three BaFei.sRuo.s As2 supercells, A, B and C 
as described in the text. 



Ru. These were with the two Ru directly on top of each 
other (A), shifted by one Fe-Fe distance (B) and shifted 
by one lattice parameter, a (C). The tetragonal lattice 
parameter for BaFe1.5Ruo.5As2 , a — 3.98 A and c = 12.95 
A were taken from Ref. [[(J and the internal coordinates 
were optimized by energy minimization. As shown in 
Fig. [21 the electronic structures for these three super- 
cells are very similar. The values at the Fermi energy for 
the three cells are N(E F )=3A7 (A), N(E F )=3A4: (B), 
and N(E F )=3.17 (C). In the following we focus on re- 
sults for supercell A, since the same conclusions would 
be drawn based on the others. 

It is found that, as the result of larger size for Ru 2+ , 
apart from expanding distances within Fe-Ru planes,^ 
the Ru-As distance (2.41 A) is also slightly larger than 
the Fe-As distances (2.30 A - 2.35 A), and accordingly 
the Ba layers are also distorted. The calculated DOS is 
shown in Fig. [3] in comparison with that of BaFe2As2- 
We did not find any additional peak associated with ex- 
tra flat bands introduced by Ru substitution. This and 
the shape of the DOS reflecting bands of mixed Fe and 
Ru character, as opposed to separate Fe and Ru derived 
bands, indicate that the electronic structure of the Fe-Ru 
system is that of a coherent alloy as might be anticipated 
from the similar electronic structures of BaRu2As2 and 
BaFe2As2- In addition, it can be clearly seen that the 
shape of total DOS near E F is almost unchanged after Ru 
substitution, except that the values decrease somewhat. 
This decrease mainly results from the enlarged d band 
width as the result of increased hopping between transi- 
tion metals and hybridization involving Ru. N(E F ) is re- 
duced to 3.17 eV _1 per BaFe1.5Ruo.5As2, as compared to 
3.28 eV _1 for BaFe2As2 (note that these GGA values are 
somewhat larger than the LDA values, e.g. N(E F )=3.06 
cV _1 for BaFe2As2 in Ref. [25h . In any case, we can con- 
clude that no additional electron carriers are introduced 
to Ru alloyed system in the sense that the band filling is 




FIG. 3: (Color online) Calculated total DOS and projec- 
tions for BaFe1.5Ruo.5As2 (top), using the quadrupled su- 
percell A (see text). The total DOS is on a per formula 
(BaFe1.5Ruo.5As2) basis while the projections are per atom. 
The total and Fe d-projected DOS for BaFe2 As2 (bottom) are 
presented for comparison. 



unchanged. 

The general shape of our density of states is simi- 
lar to that reported by Paulraj and co-workersJ. How- 
ever, there are significant differences in detail. As men- 
tioned, we find a decrease in the value of N(E F ) upon 
alloying rather than an increase. We find that for 
BaFe1.5Ruo.5As2 the distribution of Fe d states is almost 
the same as in non-spin-polarized BaFe2As2, giving the 
main contribution to the states near E F and only mod- 
estly mixing with As p states. In contrast, Ru d derived 
bands form a peak at ~-3.5 eV and overlap in energy with 
As p states concentrated within the region of -5.5 eV to 
-2.5 eV, and thus more strongly hybridized with them in 
spite of the larger Ru-As distance. This results in a re- 
duced Ru d contribution to the states near E F . Thus the 
contribution to N(E F ) from Ru d states is lower by half 
than that of the Fe d states. This plays an important role 
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FIG. 4: (Color online) Total and projected electronic DOS 
of SrFe1.5Iro.5As2, calculated in the same way as that of 
BaFe1.5Ruo.5As2. The tetragonal lattice parameters from ex- 
periment (Ref. \13tl a = 3.95 A and c = 12.22 A were used for 
supercell simulation. 

in the decreased N(Ep) for the BaFei_ a; Ru a ;As2 system 
and is a consequence of Ru-As hybridization. 



IV. SUPPRESSION OF THE SDW 

As mentioned, understanding the suppression of the 
SDW order is important. It is convenient to discuss 
the magnetic susceptibility, x(q)i which is given in the 
random-phase approximation by the enhanced Lindhard 
susceptibility x(q) = xo(q) [1 - I(q)xo(q)]~\ where 
Xo(q) is the bare susceptibility, which is a density of 
states like term and J(q) is the Stoner parameter, whose 
q dependence reflects the changing band character on 
the Fermi surface. In BaFe2As2 and other undoped Fe- 
based materials, the high N(Ep) oc Xo(0) with the large 
Stoner parameter I (~0.9 eV, characteristic of a 3d tran- 
sition element, with N(Ep) now on a per spin per Fe 
basis)^ 6 . puts them near the Stoner criterion (N(Ep)I > 
1) for itinerant ferromagnetism magnetism. 3 Further, 
the strong Fermi nesting between approximately size- 
matched electron and hole sections yields on top of this 
a peak in xo(q) a t the zone corner. This explains the 
SDW instability at the nesting vector (1/2,1/2). 

Turning to the Ru alloyed system, on the one hand, 
more hybridization with As p states and the larger size 
of the Ru Ad orbital relative to the Fe 3d orbital cause a 
reduction in N(Ep), as mentioned. On the other hand, 
the atomic-like I is reduced to the value of ~0.6 eV for 
the mid-4ci transition metal, Rui^ 6 . For this coherent al- 
loy system, the Stoner parameter is the average of the 
atomic Stoner parameters weighted by the squares of the 
partial contributions of the different d orbitals to the xo 
(i.e. N(E F ) for q=0). 26 ' 27 The rather smaller N(E F ) 



of Ru d states as mentioned would further reduce the 
contribution from Ru to the average Stoner parameter. 
Considering both / and N(Ep) it is clear that the Ru 
alloy is much less magnetic than the pure Fe compound. 
The same factors will apply to the SDW instability, and 
provide the explanation for its suppression upon Ru al- 
loying seen both experimentally, and in our direct density 
functional calculations. 



V. IRIDIUM SUBSTITUTION 

Finally, we also calculated the electronic structure for 
5d transition metal, Ir alloyed system, SrFe1.5Iro.5As2, 
As shown in Fig. [4] the substitution of Fe with Ir also 
leads to an electronic structure characteristic of a coher- 
ent alloy. Compared with BaFe1.5Ruo.5As2, the scale of 
total DOS is further reduced, reflecting still stronger hy- 
bridization and expanded band width. This is seen for 
example in the large Ir contribution to the As p derived 
bands at ~ -4 eV and the correspondingly reduced Ir 
contribution to the metal d bands from ~ -2 eV to 2 
eV relative to Ep. More importantly, Ep clearly shifts 
upwards towards the bottom of the pseudogap, which re- 
flects additional electrons that are doped into the d bands 
by alloying this Co column transition metal. The dop- 
ing amounts to one carrier per Ir since, as for Ru, the 
electronic structure is coherent and no additional bands 
are introduced. The value of N(Ep) is 2.91 eV -1 per 
formula of SrFe1.5Iro.5As2. In addition, we note that the 
Stoner parameter of Ir is rather low, I\ r ~ 0.5 eV~£ Thus, 
the suppression of SDW order in this Ir-alloyed system^ 
may be attributed to effects of both a decreased Stoner 
factor and electron doping. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we show based on density functional cal- 
culations the alloying BaFe2As2 with Ru and SrFe2As2 
with Ir results in the formation of a coherent electronic 
structure without additional localized states associated 
with the impurity atoms. As such, the suppression of 
SDW magnetic order for the Ru alloyed Ba(Sr)Fe2As2 
system is due to the increased hybridization and a re- 
duced average Stoner parameter reflecting the different 
chemistry of 3d and Ad transition elements. A simi- 
lar consideration applies to alloying with Ir, although in 
that system Ir also provides doping. We note that scat- 
tering due to Ru/Ir disorder in the actual alloys might 
also play a role in suppressing the SDW, since scatter- 
ing generally works against a nesting induced instabil- 
ity. Also it should be mentioned that the detailed role 
of Ad and 5d transition metals on the interplay between 
magnetism and superconductivity in Fe-based supercon- 
ductors might be quite complex. For example, Ru al- 
loying is found to completely suppress SDW order in 
the PrFei-^Ru^AsO system but no superconductivity 
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emerges, at least in samples investigated to date^ How- 
ever, Ir doped LaFei-^Ir^AsO is reported to show su- 
perconductivity with T c up to 10.5 Kj22 Further stud- 
ies will be helpful in understanding these differences. In 
any case, substitution of Fe with 4d and bd transition 
metals in Fe-based materials provides another avenue for 
suppressing SDW magnetic order and thereby inducing 
superconductivity in Fe-based materials, even in the ab- 
sence of doping. 
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